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Abstract 

The role played on the scalar transport by a finite, not small, correlation time, r u , 
for the noise velocity is investigated, both analytically and numerically. For small 
t u 's a mechanism leading to enhancement of transport has recently been identified 
and shown to be dominating for any type of flow. For finite non-vanishing r u 's we 
recognize the existence of a further mechanism associated with regions of anticor- 
relation of the Lagrangian advecting velocity. Depending on the extension of the 
anticorrelated regions, either an enhancement (corresponding to constructive inter- 
ference) or a depletion (corresponding to destructive interference) in the turbulent 
transport now takes place. 

PACS 47.27Qb - Turbulent diffusion; interference phenomena. 

The understanding of the mechanisms leading to transport enhancement / depletion is 
of great interest in various domains of science IJ. In the theory of scalar transport, the 
main regime of interest corresponds to the long-time behavior of the concentration field 
(e.g. the density of passive particles injected into the flow) for time scales much longer 
than the characteristic time-scale of the advecting velocity field. In this case, long-wave 
disturbances governs the slow modes of concentration and the asymptotic dynamics is 
described by a diffusive equation with renormalized diffusion coefficients, D^ , namely: 
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where (9) is the concentration field averaged locally over a volume of linear dimensions 
much larger than the typical length l Q of the advecting (incompressible) velocity field. 
The latter equation can be rigorously proved by using, for instance, asymptotic methods 
like multiscale techniques ||. From the Lagrangian viewpoint, the large-scale diffusion 
regime is a consequence of the central limit theorem, stating (roughly) that, for particle 
displacements weakly correlated, their sum tends to a Gaussian and, as an immediate 
consequence, particles undergo diffusion. 

An exact (i.e. non-perturbative) expression of the Green- Kubo type obtained by Taylor 
H allows one to express the effective diffusivities as time-integrals of Lagrangian correla- 
tions. For Lagrangian chaotic flow the correlations decay rapidly, the integral converges 
and gives a finite effective diffusivity D E . 

In many geophysical domains (for instance, to study the dispersion processes either 
in atmospheric or oceanic flows), Lagrangian trajectories, x(t), of particles advected by a 
prescribed velocity field, U(x(t)), and plunged in a noise, small-scale velocity, u(t), are 
properly modelled by the Langevin equation: 



where the small-scale turbulence, u(t), is described by a Markovian process with finite 
correlation time, r u , rather than by a white-noise process, on account of experimental 
observations |4], |5j suggesting the presence of some form of finite correlation for the smallest 
scales of motion. 

In the following, the advecting velocity field U(x(t)) will be assumed incompressible, the 
noise term u(t) a Gaussian random process having zero average value and the colored- 
noise correlation function given by: 



where D e is the isotropic eddy-diffusivity arising from the smallest scales of turbulent 
motion (not explicitly solved) and r u is their correlation time. 

The main question addressed in the present Letter, concerns the role played by the 
finite (not necessarily small) correlation time, t u , of the noise velocity in the large-scale 
dynamics and thus in the behavior of the effective diffusivity, D E . Specifically, we are 
interested in the difference 



for small, fixed, D e and different values of t u . The argument of D E is meant to stress 
that the eddy-diffusivity depends on both D e and r u due to the dependence of Lagrangian 
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trajectories on the colored-noise, as a whole. 

Physically, positive (negative) values of I correspond to constructive (destructive) inter- 
ference induced by the dependence of Lagrangian trajectories on the colored-noise. 
In Refs. || and j7| the interference problem was investigated focusing the attention on 
the variation of the sign of I for small D e and r u = (i.e. for white- in-time processes). 
It was conjectured in Ref. || that the interference process should always be destructive, 
and actually shown in Ref. [7] that it can be either constructive or destructive, depending 
on the extension of the anticorrelated regions of the velocity field. 

The effect of a small t u on the variation of I was investigated in Ref. M. Accordingly, 
splitting I into two distinct parts, / = Ipe + I Tu , the first one being the interference 
contribution of D e when t u = 0, and the second the correction due to a non- vanishing 
t u , it was shown perturbatively in the parameter t u /t v <C 1, that I Tu > for any type 
of flow, Ty being the correlation time of the advecting component U. Physically, the 
mechanism leading to positive interference, and thus to enhanced diffusion, was identified 
in the augmentation of coherence in the diffusion process induced by the small correlation 
time t u . Indeed, the latter makes the particles of the diffusing substance to forget less 
rapidly their past than in the presence of white-noise turbulence. This fact increases the 
Lagrangian correlation time and therefore the eddy-diffusivity. 

In this Letter, we want to remark that the above mechanism is certainly present for all 
values of r u , but when the latter is large enough another effect can be even stronger. To 
show this, let us consider the case when the Lagrangian correlation function presents neg- 
ative correlated regions. The contribution given by such regions in the integral along the 
trajectories required in order to obtain the eddy-diffusivities is negative and a reduction 
of the diffusion may occur. The role played by t u is thus twofold. On one hand, according 
to the results of Ref. ||, it permits the particles to go away from their initial positions 
more slowly. This effect increases the Lagrangian correlation time and therefore D . On 
the other hand, it works to maintain the particles inside the trapping zones where the 
velocity is anticorrelated. The slowing-down due to these anticorrelated regions is thus 
enhanced and diffusion is reduced. Due to the well-known sensitivity of transport to slow 
parts of the trajectories, the latter mechanism can become stronger than the former and 
a reduction in the transport may occur (i.e. I Tu < ). 

In the following we show that the destructive effect associated with anticorrelated 
regions can be stronger than the constructive effects. To that end, in order to simplify 
the problem, let us consider a very simple model for the advecting velocity, U, whose 
stream function is defined as: 

^{y, t) = 4>o [a(t) cos (k y) + (3(t) sin (k y)] . (5) 

For this 2-d random parallel flow, as in Ref. J7j, the (random) processes a(t) and /3(f) are 
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assumed Gaussian, independent and with auto-correlation function given by 

C(t) = e- |t|/ru cos (ut) (6) 

where t v is the correlation time of U and u is the parameter controlling the extension of 
anticorrelated regions. The form of @ permits us to obtain a finite and non-vanishing 
effective diffusivity also for D e = and, moreover, represents the simplest way to introduce 
recirculating zones where the scalar tends to be homogenized and trapped. 
The following example of cellular flow defined as 

U(x,y) = U (cos(2/),cos(:e)) , (7) 

the well-known ! BC flow ||, [IT], |TTJ whose streamlines form a closed structure made of 
four cells in each periodicity box, shows indeed that auto-correlation functions like (||) are 
associated to recirculation. In order to show this point, we have integrated the Langevin 
equation by using a second-order Runge-Kutta scheme, with the advecting velocity field 
given by eq. (0) superimposed to a white-in-time noise. The averages necessary to evaluate 
the auto-correlation function C{t) = (U(x(t)) ■ U(x(0))) /U 2 , (U 2 = (U-U) is the energy 
of the advecting velocity field) have been made over different realizations and performed by 
uniformly distributing 10 6 particles in the four cells of the periodicity box. The obtained 
auto-correlation function is shown in Fig. 1. The qualitative resemblance between the BC 
auto-correlation function and that given by the simple expression makes meaningful 
the use of the latter to mimic the statistical effect of recirculation in the diffusion process. 
Notice that, although possible, the main disadvantage in studying the interference problem 
by using the BC flow is that only numerical analysis seems to be feasible. This is not the 
case for the random parallel flow (|5p with the auto-correlation (|(J), where the interference 
mechanisms can be investigated analytically. To that end, we use asymptotic methods 
since we are interested in the dynamics of the field 9 on scales of the order of L, larger 
than l , the typical scale of the advecting velocity field, U. The ratio l /L ~ 0(e) 
with e C 1, the parameter controlling the scale separation, naturally suggests to look 
for a perturbative approach. The perturbation is however singular El and multiscale 



techniques [|] are thus needed to treat the singularities. For the random parallel flow, the 
auxiliary equation (a partial differential equation which needs to be solved to obtain the 
eddy-diffusivity) can be tackled analytically and a simple expression for the component 
of the effective diffusivity tensor parallel to the direction of the velocity U (the normal 



component being equal to D e ) can be obtained |L3[ and its expression reads: 
D E (D e ,T u )=D e + u 2 j C(t) e * c 



t+T u [e 



dt 



where Uq = ipo^o and C(t) has a typical correlation time t v as in 
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Figure 1: The velocity auto-correlation function (U(x(t)) ■ U(x(0))) /U 2 as a function 
of the time t for the 'BC flow ([?]). The curve has been obtained by integrating, with 
a second-order Runge-Kutta scheme, the Langevin equation for the coordinate of scalar 
particles in the presence of a white-in-time noise. 



The two opposite limits t u /t v <C 1 and t u /t v 1 can be easily investigated from (§). 
About the former, it is immediately checked that the results obtained in Ref. || are here 
reproduced by expanding at the first order in both t u and D e the term inside the curly 
brackets in the r.h.s. of (§). As a result, the interference term I Tu reads: 



D e k'uir u \ C(t) dt . 
Jo 



(9) 



The second limit can be easily exploited by expanding the same term inside the curly 
brackets but now for t u /t v ^> 1. The expression for D E corrected at the order 1/t u is 
thus obtained: 



/ C(t) 
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D E = D e + ul 

from which the interference term / immediately follows: 

ulk 2 D e r°° 2 „. . , 
I = — 5_5 — / t 2 C(t) dt . 
2t u Jo 



dt 



(10) 
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It is easy to verify from (|8J) that, for both t u = and small D e (i.e. D e k r ^ j r v 1), the 
interference term associated to D e is given by: 



L£)e 
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and thus, from (0), it follows: 



(1 poo poo \ 

— / t 2 C(t)dt- tC(t)dt) . (13) 
2t u Jo Jo J 

For large r u 's, the character of the interference due to the effect of t u on the Lagrangian 
trajectories is thus related to the variation of I Tu with r u . 

By inserting the expression (|6|) for C(t) into (p~3[ ) we obtain, after simple integrations, the 
following condition for destructive interference: 

I Tu < when ur v > 1 . (14) 

According to our previous discussions, destructive interference appears when the extension 
of anticorrelated regions (controlled by uj) is large enough. Such phenomenon, found in 
the limit of large r u 's, is actually present for a wide range of values of r u 's provided that 
T u /r v <jt 1. 

To show this point, we have computed numerically the integral in the r.h.s. of @ from 
which we have subtracted the contribution Joe + D E (0, 0) to obtain I Tu . The behavior of 
I Tu as a function of ujt v for t u /t v = 0.05 (long-dashed line), 0.5 (dashed), 1.0 (dotted) and 
5.0 (dot-dashed) is shown in Fig. 2 for k\D e r^ = 0.05 and Uq = 1. From this figure we 
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Figure 2: The interference term I Tu as a function of ujt v , for different values of t u /t v : 0.05 
(long-dashed line), 0.5 (dashed line), 1.0 (dotted line) and 5.0 (dot-dashed). Curves have 
been obtained by numerical computation of integral (ffl) with — 1 and \^D e r v = 0.05. 

observe that the smaller t u /t v , the larger the values of ut v at which the crossover from 
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constructive to destructive interference takes place. For all r u 's, we can then conclude 
that the condition (|TJ|) should be replaced by: 

I Tu < when ujt v > K c (15) 

with K c depending on t u /t v , K c > 1 (on account of (|14|)) and lim Tti ^ K c = oo (on account 
of the positive character of the interference when t u — » 0). Condition ( pT5| ) is corroborated 
by the results of Fig. 3, where the values, K c 's, of ujt v at which the crossover takes place, 
have been plotted for different values of t u /t v . 
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Figure 3: Values K c of uot v corresponding to the crossover from constructive to destructive 
interference (i.e. from I Tu > to I Tu < 0), as a function of the ratio t u /t v . Values of K c 
reported in figure have been obtained by numerical computation of integral ([8]) with the 
same parameters of Fig. 2. 

In conclusion, the effect on the scalar transport due to a finite correlation time t u of 
the noise velocity has been investigated. We have shown, both analytically and numeri- 
cally, that for t u large enough two mechanisms exist and act in competition to determine 
the interference character. The first mechanism works to enhance the transport due to 
an augmentation of coherence induced in the diffusion process. The second mechanism 
maintains particles inside the zones of trapping (where the auto-correlation of the advert- 
ing velocity field is negative) for a time longer than the time relative to the white-noise 
case. Depending on the extension of anti-correlated regions, constructive or destructive 
interference take place. For finite values of r u , the symmetry with the scenario presented 
in Ref. is thus perfectly restored. 




7 



Acknowledgements We are particularly grateful to M. Vergassola and A. Vulpiani 
for very stimulating discussions and suggestions. Very useful discussions with A. Crisanti, 
R. Pasmanter and E. Zambianchi are also acknowledged. The work of AM has been 
supported by the "Henri Poincare" fellowship (Centre National de la Recherche Scien- 
tifique and Conseil General des Alpes Maritimes). PC is grateful to the European Science 
Foundation for the TAO exchange grant. PC has been partially supported by INFM 
(PRA-TURBO) and by MURST (no. 9702265437). 

References 

[1] H.K. Moffatt, Rep. Prog. Phys., 46, 621 (1983). 

[2] A. Bensoussan, J.-L. Lions and G. Papanicolaou, Asymptotic Analysis for Periodic 
Structures (North-Holland, Amsterdam, 1978) 

[3] G.I. Taylor, Proc. Lond. Math. Soc. Ser. 2, 20, 196 (1921). 

[4] R.E Davis, J. Geophys. Res., 90, 4756 (1985). 

[5] P.M. Poulain and P.P. Niiler, J. Phys. Oceanogr., 19, 1588 (1989). 

[6] PG. Saffman, J. Fluid Mech., 8, 273 (1960). 

[7] A. Mazzino and M. Vergassola, Europhys. Lett., 37, 535 (1997). 

[8] P. Castiglione and A. Mazzino, Europhys. Lett., 43, 522 (1998). 

[9] VI. Arnold, C.R. Acad. Sci. Pans A, 261, 17 (1965). 

[10] M. Henon, C. R. Acad. Sci. Pans A, 262, 312 (1966). 

[11] T. Dombre, U. Frisch, J. M. Greene, M. Henon, A. Mehr and A.M. Soward, J. Fluid 
Mech., 167, 353 (1986). 

[12] CM. Bender and S.A. Orszag, Advanced Mathematical Methods for Scientists and 
Engineers (McGraw-Hill, 1978). 

[13] P. Castiglione and A. Crisanti, "Dispersion of passive tracers for a Langevin equation 
with non delta-correlated noise", submitted to Phys. Rev. E, (1998). 



8 



